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Abstract 

It’s a common knowledge that the quasielastic neutrino-neutron 
and antineutrino-proton cross sections tend to the same constant as 
(anti)neutrino energy becomes high. In this paper we calculate the 
exact expression of the limit in terms of the parameters describing 
quasielastic scattering. We check that even at very high energies only 
small absolute values of the four-momentum transfer contribute to 
the cross section, hence the Fermi theory can be applied. The dipole 
approximation of the form factors allows to perform analytic calcula¬ 
tions. Obtained results are neutrino-flavour independent. 

1 Introduction 

Quasielastic neutrino scattering plays a dominant role in neutrino-nucleon 
reactions at energies below 1 GeV. When neutrino energy increases another 
channels open and quasielastic processes become less important. At high 
energy the total cross section for neutrino scattering is approximately pro¬ 
portional to the value of the energy while the quasielastic cross section is 
roughly constant. The latter behaviour is known on the basis of numerical 
computations but as far as we know, it has not been shown analytically yet. 

The quasielastic cross section is usually calculated within the Fermi the¬ 
ory. At low energies four-momentum transfer is understood to fulhl the 
condition |g^| -C where Mw = 80.4 GeV is hF-boson mass. It will be 
shown in Sec. 01 that in fact even for very-high-energy neutrinos overwhelm¬ 
ing contribution to the cross section satisfies such constraint, therefore the 
use of the Fermi theory is well justified. 
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Radiative corrections are not taken into acconnt, but in Sec.Elwe estimate 
that they can be neglected. 

In the theoretical description of the neutrino-nucleon interaction the ha¬ 
dronic current is expressed in terms of the four form factors due to Lorentz 
invariance and assumption that there are no second-class currents. The form 
factors can be expressed in various ways, see jH]. We consider dipole form 
factors because of their simplicity in analytic calculations. 

The quasielastic cross section for neutrino-neutron scattering can be writ¬ 
ten as P 


a = 


M^G% cos2 Oc 


dq^ 


A{q^) - Biq^) 




M2 


+ W)- 


M4 


( 1 ) 


where M 

A{q^) = 


= (m„ -|- nip)/2 is the average nucleon mass and 



- ^ +nw + + mFAFp) + jfi\Fpf) 


B(f) = -AmFl, + iF^)FX), 




In above formulae mi is charged-lepton mass, — neutrino energy and 
i = qip — iin — where [ip and pn are the proton and neutron magnetic 
moments respectively. In the case of antineutrino-proton scattering —B{q^) 
in Eq. (HD should be replaced by +B{q^). We also need to know the interval 
of integration [{q^) a, {q^) b\- 


{q^)A = 

iq% = 


mfiE^ + M) - 2MEl - ^ 
2E^ + M ’ 

mfM 

m‘f{E^ + M) - 2MEI - 


( 2 ) 


with Z\ = (2MEl - m/E^)^ - Am/M^E/. 

As it was mentioned before, in this paper we will consider dipole form 
factors. Using the Sachs form factors 


GW) = 


(l-g2/M2)2’ 


GW) = 


1 + e 


(l-g2/M2)2’ 


2 

















Gp 

1.1803 

10-VGeV2 

cos 9c 

0.9740 


9a 

-1.267 



3.7059 

fJ-N 

Ma 

1.001 

GeV 

M2 

0.71 

GeV2 


Table 1: The values of the constants used in numerical calculations 
the vector form factors can be expressed in the following way: 


(f 






4M2 

Vt2\ , 


whereas the pseudoscalar form factor Fp is related to the axial one due to 
PCAC hypothesis: 


FA{q^) 


9A 

(l-g2/M2)2’ 


Fp{q^) 


2M^FA{q^) 
ml — g2 


By we denoted the pion mass. 


2 High-energy limit 

If neutrino energy Ei, is high enough to fulhl the condition /Ey <C 1, 

where = ma.x{mi, M, My, Ma}, one can write 


Z\ = AM^E'^ 


p2 rnt 


„ mj 
Ey + - m. 


what results in 


{q^)A 


-2MEy, 


AM^E: 


2 7714 


{q^)B - 0 . 

The cross section m is the sum of terms 

6^ 


(3) 


a = 


-/? = 


K = 


AEl 

Q 

AEl 

Q 

AEl 


dq^M^A{q^), 
dq^B{q^) (s -u), 


( 4 ) 


dq^C{q^ 


u] 


M2 
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Figure 1: The cross sections’ dependence on neutrino energy, aoo stands 
for the high-energy limit of a calculated in this paper. Experimental data 
for quasielastic scattering from D 2 target are taken from ANL 1973 [3], 
ANL 1977 0, BNL 1981 0, FNAL 1983 0 and CERN-WA25 1990 0. 


where we have introduced the compact notation for the constant factor 

Gl cos^ Oc 
^ ~ 2ti 

The hrst term, that is a, tends to zero as neutrino energy becomes inhnite. 
We will show it in Appendix ^ 

Next, in Appendix |Bl it is calculated directly that in the discussed limit f3 
also approaches zero. This result is consistent with Pomeranchuk’s theorem 
(see generalization in 0 ), which states that as Ej, —>■ cx) the neutrino and 
antineutrino cross sections become equal. They differ in sign of the (3 term, 
so the term should tend to zero. 

Thus only k, gives a nonzero contribution to the high-energy (anti)neu- 
trino quasielastic cross section: 


(Too = lim a = lim k. 

Ei/^OO Eu^OO 


Our main result can be can be written in the form 
cos^ 6c ■ 


o-o 


Qtt 


Ml- + g\Ml + - Ml) 


(4M2 - m2 )2 
3e(e + 2)M8/ AM- 
(4M2 - m2 )3 Um2 - 


’ 4M2 

K W 
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Figure 2: The dependence of the high-energy limit of the cross section cToo 
on the axial mass. Marked point represents the value of as in jH]. 


Detailed calculations are presented in Appendix O Introducing notation p = 
AM'^/My and /i = ^ -|- 1 we can write it in the more compact way: 


(To 


G'^p cos^ 9c 
Gtt 


+ g\Ml + 2M^ 


3M 


— 1 

2 - 1 
V 


1 - 


P 


(p- l)Hp- 1 


P 


ln(p) - 1 


(5) 


The above expression does not depend on the charged-lepton mass, therefore 
the Ey ^ oo limit of the cross section is equal for all the neutrinos and 
antineutrinos, see also Fig. [U 


3 Discussion 

To obtain numerical value of the limit we assume values of the constants for 
dipole form factors as in [Hj, see Tab.Q Note the corrected value of the axial 
mass: Ma = 1.001 ± 0.020 GeV. Then 

Uoo = 0.956 X 10"^® cml 

We observe next that none of the four terms in Eq. can be neglected. 
Contribution of the term with the axial form factor is equal to about 46%. 
The dependence of cxoo on the value of the axial mass is shown in Fig. El 
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Figure 3: The ratio of the difference between the cross section aw with 
the hF-boson propagator and the Fermi-theory cross section normalized with 
respect to aw itself. 

It is necessary to check if our approach based on the Fermi theory is 
consistent. We do it numerically by computing the cross section with the W- 
boson propagator aw and comparing the result with the cross section within 
the Fermi theory a. Fig. IHl presents the dependence of the ratio 

^ _ o'w — ^ 
aw 

on neutrino energy. When E^y > 50 GeV the ratio R is roughly constant and 
less than 0.01% (for each flavour). It means that only small four-momentum 
transfers \q^\ contribute to the quasielastic cross section, thus calculations 
within the Fermi theory are reasonable even for very high neutrino energies. 

In calculations of the limit of the cross section no radiative corrections 
were taken into account. We guess that corrections to quasielastic scattering 
are of the same order of magnitude as to deep inelastic scattering, i.e. they are 
roughly constant and of the order of half a percent [H] (the value refers to the 
corrections which comes from bremsstrahlung of the charged lepton, W boson 
and quarks). If the hypothesis is true, it makes them of low importance unless 
experiments reach very high precision. 

More important improvements could come from the non-dipole form fac¬ 
tors as in [H]. Presented there hgures suggest that they would yield the value 
of the limit 3% smaller with respect to our result, but unfortunately “BBA- 
2003 Form Factors” are practically unapplicable to analytic calculations. 
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A Why a tends to zero 

The a term defined in Eq. m is an integral of rational fnnction of divided 
by nentrino energy sqnared. As ^ oo, a wonldn’t tend to zero only if the 
integral rose at least as E^. The form of the limits © implies that the lower 
one always gives zero and only the npper one could produce nonzero terms, 
if the integrand is of the order at least one in Let’s write explicitly the 
term of the highest order for each form factor, keeping in memory that Ep 
can be expressed by Ea- 

glM\ {q^f 


4VM2y 

A\M^) 


\EA^ = 


ipl |2 _ 

^v\ - 


4M4 {Ml-q 
2 


2'|4’ 




{q^Y (4M2-g2(^ + l))^ 




4M4 (4]pf2 _ _ q2^ ‘ 

eM^ 1 (g')^ 


i2'|4 


M2 - q\ 

We can see that each one of them is a proper fraction, so as neutrino energy 
becomes inhnite a tends to zero. To be sure, let’s perform the calculation 
for the second of above expressions. We can obtain easy-to-integrate form 
by decomposing it into partial fractions: 


h—V 

4 VM2/ 


I pi |2 _ 

Tvl - 


MS r 
4M4 


+ 0{l/q^ 


.{AM^-q^) (M^-q^) 

where c is a constant and 0(1/q^ denotes terms of lower order in g^. As 
neutrino energy becomes high the limits of integration are given by © hence 


1 

4 




^l|2 

yl 


MS r 
4M4 


, My , 

2c In + other constants 


The above integral tends to a constant as E^ ^ oo. Only higher order term 
in g^ could give result increasing with E^y but there is no such term in a. 
Since that for the whole expression holds true that 

const 


a 


E?, 


0 . 


B Why f3 tends to zero 

In the frame in which target nucleon is at rest {s — u) = {AME^ — rrif + g^), 
so the quantity dehned in Eq. 0 can be explicitly written as 

ligAG{MAMvY ^ 




AM'^El 


dq^(^‘B{AME,y — rnf) + ®g^). 
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where 


« = 


T 


To perform the integration one need to decompose the integrand into partial 
fractions: 

Ml + Ml 


® = 


= 


7?2 




Ml 


+ 


Ml 


(M2-g2)2 


+ 


Ml 


+ 


Ml 


(Ml-g2)2 [Ml-q^f 


+ 


TZa 

2MlMl 

TZa 


Ml - g2 Ml- g2 
1 1 


Ml-g2 


Ml 




where TZa = Ml — Ml. As -C 1, after integrating in the limits (jSI) 

we obtain 


^ ^ 2figAg{MAMvY ^ Ml + Ml ^^My 


ME,1Z\ 


TZ, 


Ma 


0 . 


C Why K, tends to constant 

The last term in Eq. expressed by the form factors is 


K = 


{AME^ 


dq^ 


^l|2 

y| 


q 


4M2 


\iEll + \EA[ 


U] 


For convenience we separate the axial part from the vector one: 

^ f 1 2 ITT \2/^ 


Ka = 




{AME^ 


{AME^ 


dq^\EAl{s-uy, 
,2 


dg2 


I pi |2 _ q 

I 4M2 


m 


2|2 


s - u 


To evalnate the integral 


9aGMI /• , 2 (g-^) 

^ (4ME,)2 J ^ {Ml- q 


214' 


one needs to know decomposition of the integrand. If we add Ma and —Ma 
to (s — u) and sqnare it in the following way 

(s-n)2 = {AME,-m^ + Miy-2{AME,-ml + Ml){Ml-qy + {Ml-qyi 


we will get 


u 


{AME^ - ml + Mil 2{4ME^ - ml + Ml) 


{Ml - q^Y 


{Ml - qY' 


(Ml - «2)3 


+ 


{Ml-q^f 



































It means that as neutrino energy fulfils condition <C 1, integration 

in the limits © leads to 


Ka 


g\QM\ ( 2mf + Ml 


V 


and 


lim Ka = Q 

Ey^oo 3 


AMEy 

9\MI 


The integrand in dehnition of Kv, he. 


IF, 


1 |2 
yl 




4M2 


mv = 1 




2 N -1 


4M2 


2 N -4r 




M^J 


h'(l- 


4M2 


+ 1-F 


with /r = .^ + 1, can be written as 


IF 


1 |2 
ul 


9^ |.p2|2 

14-hyl = 


4M2M8 (/i2 - 1) 


4M2 


(4M2-g2)(iVf2 -g2)4- 


Let’s denote the last-fraction’s numerator as /C = 4M^My(yU^ ~ !)• Above 
expression decomposed into partial fractions is 


I f 1|2 _ _^_|yp2|2 _ a: 


4M2' 


n\,\Ml-q^ AM^-q^J nl{Ml - q^f 

K, IC„ 


+ 


F 


nUMl-qY^ nv{Ml-qY^^ 


where TZy = AM‘^ — My and = My{AM‘^ — iYM^). By repeating the 
trick made during the computation of Ka we obtain 


T7i |Z 

^v\ - 


AM^ 


mi {s-uY = 


Cl 


Cl 


+ 


C2 


M^-q^ AM^-q"^ [M^ - q^Y 

C 3 , C 4 


+ 


+ 


(M2-g2)3 (M2-g2)4’ 


where coefficients are: 

F 

n 

■'V 

) 

C2 = ^ 

1 


c, = mAMEy-m^ + AMY^ 

Ky 

F 2 


nl 


2\2 


F 


- (4MF,, - mi + AMY 


C 3 = 


nl 

F, 


K,(AMEy - mi + Ml 


2 

F 


(F^F 


yj 


12 -, 


F 


C 4 = ^(4MF,.-m2 + M2)2. 
/vy 
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For neutrino energy Ey 3> we conclude that integration over {dq^) 

leads to 


Q 

{^MEyY 

<3 


2 / *^1 


)2 j ^ 4M^-qy 


GE , 4M2 4M2 - m? 




g / , 2 ^ g/C / 4M^ - 

(4ME,)2 7 - g')' ^y^y V ^ 2MEy 


;, 2 C3 g ( E{M^-m^)-E,TZv 

^ (M2 - g2)3 ^ 2M47^^ V 2ME^ 

' , 2 C4 GEf, ( M^-mf\ 


{AMEy 


^ / j 2 '-4 _^ , ^’--‘y 

{AMEyf J ^ (M2 - g2)4 3M6 7^y V 2ME, 


The K term is the sum of Ka and k^, therefore 


,. 4M2 3/C 3/C 

e!^oo^ ^ 3 ^ 37^y iTZ^ M2 M2 7^^ ^ 2M4 7^y ^ M^.' 


Recall that /C^ = My{AM‘^ — j^My) and IZy = 4M2 — My, hence we obtain 


^ = M7(4M2 - /i2M7) = M77^y - {^y - 1)M7. 


Next, constant factor /C = AM'^My^f^ — 1), so 


2(M2 - m2 ) + 3M2 


I y.2 


It means that the limit of the cross section is equal to 


lim a = 

Eu^OO 


GlcoG^dc\..2 , 2 


[^v + 9aMa+ [M 


^ (4M2-M2)2' 

3(/i2 - 1)M^ / 4M2 


^ _In _ i 

{AM^ - M^yXAM^ - M2 J\' 


Denoting p = AM'^/My we can write this formula in the following way: 


lim a = 

Eu^OO 


Gp cos2 9c 
Gtt 


A4 + j>q + 2Af"(^^(i-i) 

, 0..2 / P 
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